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Abstract: 

We compute the effects induced by the use of small CMB maps on the measurement of 
the Ci coefficients of the angular power spectrum and show that small systematic effects 
have to be taken into account. We also compute numerically the cosmic variance and 
covariance of the Ci spectrum for various spherical cap like maps. Comparisons with 
simulations are presented. The calculations are done using the standard method based on 
the spherical harmonic transform or using the temperature angular correlation spectrum. 
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I Introduction 



The field of cosmic microwave background (CMB) anisotropics has dramatically ad- 
vanced over the last decade especially on its observational front. Satellite experiments 
(COBE,WMAP) observed the whole sky. But to study the high multipoles of the power 
spectrum, balloon-borne experiments have emerged. 

The BOOMERANG ( [De Bernardis et al, 2000] ), the MAXIMA QHanany et al, 2000| ) 
and the ARCHEOPS ( |Benoit et al, 2003| ) balloon experiments made measurements of 
the first doppler peak of the CMB spectrum. Due to a limited observing time and technical 
constraints, these experiments only observed a small fraction of the sky. BOOMERANG 
observed a region of 1800 deg 2 and measured the spectrum up to multipoles l max ~ 1000 
and MAXIMA observed 124 deg 2 and measured the spectrum up to l max ~ 785. In the 
future, the OLIMPO experiment (|M asi et al, 2006] ) will observe 300 deg 2 and measure 
the spectrum up to l max ~ 2500. 

In this paper we compute the effects induced by the use of small CMB maps on the 
measurement of the C\ spectrum at high multipole moments. Two methods are used to 
compute the C\ coefficients. The first one uses the Fast Fourier Transform (FFT) technics 
and the second one the angular correlation spectrum. These two technics are complemen- 
tary and lead to different biases. We show how the variance behaves with respect to 
the map size and how the apodization works in the case of the angular correlation spec- 
trum. These calculations will be compared with simulated CMB maps of various sizes 
and shapes. Numerical calculations provide a fast tool to assess most of the biases. 
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II Estimating the power spectrum: formalism 
for the full sphere. 

ii. 1 Spectrum and angular correlation fonction 
definition. 

The observed temperature T on the 2-sphere is a random field that can be expanded 
on the spherical harmonic basis: 

oo +i 



t{9A) = E E aim^rCM (n.i) 

l=o m=-l 

where the ai m coefficients are random variables. With f2 = (9, 0) one haa3 

a lm = J T(n) r™(n) dn (11.2) 

T(Q) is a real field so ai- m = {— l) m ««m which is a consequence of the relation 

Y~ m ^ -^mym 

The observed sky is a particular realisation of that random field. If one assumes uncorre- 
lated ai m coefficients and isotropy, one has: 

< a lmO'l'm' > ens u l I' u m m 

Ci (II.3) 

where the Ci are specified by the cosmological theory and the symbol <> ens means aver- 
aging over many sky realisations. 

The temperature angular correlation fonction is defined as: 

n') = < t(q)t(q') > ens (ii.4) 

where Q , Q' are the directions of measurement. The temperature angular correlation 
function is related to the spectrum of the primordial fluctuations. With the isotropy 
hypothesis, the angular correlation function is only a function of the angular separation 

7 with cos(7) = Q .Q' . We have: 

£(7) = < T(n) T(ft') > ens 

= EE <a^a l , m/ > ens Y^(n)Y l 7 l '{n') 

l,m l',m' 

= E E ci Mw W^W^') 

l,m l',m' 
l,m 

= £ E(^ + 1) Q P,(cos( T )) 
I 

This relation can be inverted and finally: 



The overbar means complex conjugation 
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= ^E( 2/ + 1 )^ p ^ cos W) 

(II.5) 



47T 

1=0 



Q = 2tt J £(7) i^(cos(7)) dcos( 7 ) 

7=0 



In the rest of these notes, the isotropy hypothesis (III. 3D will always be assumed. 

ii.2 Power spectrum estimator. 

In practice, one observes a unique sky realisation and estimators of the power spectrum 
and the angular correlation function have to be constructed: instead of averaging over 
many sky realisations, we can use the sky isotropy hypothesis. The bias and variance of 
such estimators have to be computed. 

We can have an estimation of the ai m by using ( 111.21) . Then we get an estimator C\ of the 
Ci spectrum by averaging over m: 

i 1 

Q = Yl I a ^ I' ( IL6 ) 

m=— I 

For the whole sky, this estimator is unbiased: 

j +i I +l 

<c l>ens = — x <m m -a lm > ens = — - J2q = a 

m=—l m=—l 

We can compute the variance of the estimator: 

V(Ci,Ci>) = < C\Cy > ens ~ < C\ >< Cy > ens 

(2£+l)(2Z'+l) ^ O'lm^lm^V 'm' '0>V 'in' ^ > ens G\Gy 

mm' 

Provided that the a\ m are gaussian random variables, their fourth order moments can be 
expressed with their second order moments and we get: 

V{Ci,Cy) = (2 ; +1 w 2; / +1 ) J2 {ClCy(l + 5u>5 mm >) + ( — l) m + m CiCy5u>5 m - m >) — CiCy 



mm' 



(2«+l)(2/'+l) 
1 

(2Z+1)(2Z'+1) 

mm 

i or '2 

(2Z+1)(2Z'+1) 

mm 

2m 



(ClCy + C l 8iy8 mm i) + CfSu'S 

m—m ) - QCy 

nm' 

■mm' 
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Thus the Ci are unbiased and independants, they follow a x 2 distribution with 21 + 1 

degrees of freedom and their variances are: 

°h - si < IL7 > 



Figure [T] shows an example of a power spectrum C\ and the variance of its estimator Ci 
obtain with the CMBFAST code QSeljak et al, 1996| and |Zaldarriaga et al, 2000| ) with 



the standard set of cosmological parameters. 

ii.3 Angular correlation function estimator. 

In order to obtain an estimator £ of the angular correlation function, we can average 
over all pairs of directions Q , Q' of the observed sky (keeping the angle 7 between them 
fixed) : 

= J T(fii)T(n 2 ) dn 1( in 2 6(n[.n 2 - cos( 7 )) 

S 2 xS 2 

with 

A%) = J d^xd^a 5<$h£h ~ cos(7)) 

S 2 xS 2 

As T(fi) = E a lm Yr(tt) , we get: 

l,m 

1(7) = J-t / ^T^O^M dOidOa - cos( 7 )) 

' h,mxl2,rri2 S 2 xS 2 

The internal integral X is computed in appendix |A] and leads to the estimator value: 
£W = a/TyT ^ ah^ra« 2m2 2vr 5 il/2 <5 mim2 P i2 (cos(7)) 

= ATM; 

Averaging over realisations gives us: 

<i(l)>ens = <Ohn~ai m > ens 2?T P(cOs(7)) 

1 

AT( 7 ) 



7TZ)J2 aim 27T P;(cOs( 7 )) 



l,m 

7 J,27rE(2/ + l)C,fl(oos( 7 )) 



Now A/"( 7 ) has to be computed. A/*( 7 ) is equal to the previously computed integral with 
T(fi) = 1. So, replacing T(fi) by v^47r^o°(^)' or = V^^zo^mo and C7 = A-k5iq. As 
Po(cos(7)) = 1, one finally obtains A/"( 7 ) = Af = 8n 2 and: 



< >ens = £( 7 ) = E ( 2/ + ^(C0S( 7 )) (II.8) 



47T 

The estimator £(7) is unbiased. 



Figure [T] shows the angular correlation function £ (7) and the variance of its estimator 
1(7)- 

Defining # = (^) , the covariance of the estimator is: 

£ (< >ens - < |(7) >ens< >ens) 

= EE -^(7) Pl'il') < a lmWm a l'm'Wm' > ens 
Iml'm' 

— EE Pl\l) Pl'\l) < a lmdlm >ens< dl'm'Wrn/ >ens 
Im I'm' 

As the temperature field is real, and provided that the a\ m are gaussian random variables, 
one has: 

£ (< |(7)I(V) >ens ~ < £(7) >ens< £(Y) >« 

= EE ^(7) W)acv(i + <^mmO + EE ^(7) M-fiCifcSr. 

Iml'm' Iml'm' 

-EE ^(7) WW 

Imi'm' 

= 2 EC? fl( 7 ) A(V) 
Finally, by replacing by its value, we obtain: 



< 1(7)1(7') >ens - < |(7) >ens< £(Y) >en S = (^E( 2/ + ^ ^W) (H-9) 

The estimators of the angular correlation function values 
in different directions are correlated. 



Figure [T] shows the correlation matrix of the angular correlation estimator £(7). 

ii.4 What about computing the C\ on a portion 
of sphere ? 

There are two ways to compute the C\ coefficients. Given a temperature field on 
the sphere, one can compute the ai m with (111.20 and then the C\ using the unbiased 

estimator fllL6|) . Equation AO) is: a im ~ / P z m (cos(^))rfcos(0) f T(0, 0)e" im ^0. 
The inner integral is a Fourier transform which can be fastly computed with the Fast 
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Fourier Transform (FFT) technics (see Natoli et al, 1997] ). Note that this is possible only 
if the pixels lie on isolatitude lines. The other way is to compute the angular correlation 
function £(7), which is unbiased, and then use the second formula in (111.5p . 

If the observed map is not the full sphere we can adapt these methods. For the 
first one we can still compute pseudo a>i m by setting the temperature to zero outside the 
observation zone. This leads to a systematic bias that will be studied in section IIHI 
This also introduces correlations among the measured Ci and increases their variances. 
Numerical calculations can be performed in the case where the map is a spherical cap. 
This is explored in section II III and its last section presents detailed comparisons with 
simulations. This will show that the variance does not depend very much on the map 
shape. The evolution of the variance with the map size is interpreted in this section. 

In section HVl we shall discuss the use of the angular correlation function £(7) measured 
on the observed map (so the angular separation 7 is limited to Ju m ). We will show that this 
quantity is still unbiased, as one intuitively guesses. One can still use the formula ( 111.51) to 
extract the C\ but, since the correlation function is defined only on some interval [0, ju m ], 
the Legendre polynomials are no longer orthogonal. This introduces wild fluctuations in 
the C[. This is cured by "apodizing" the angular correlation function and will be described 
in section IV! 

Four appendices explain technical details to make these notes self contained. 
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Figure 1: 

Top left: power spectra versus the multipole moment 

Top right: angular correlation function versus the separation angle (deg) 

The red curves shows the dispersion of the estimator due to cosmic variance (values have 

been mutiply by 10 for readability) 

9¥ Bottom: the correlation matrix of the angular correlation function estimator. 
The upper left corner is (0, 0) degree and lower right corner is (180, 180) degres. 
The LUT ranges from —1 (dark) to 1 (white). 
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Ill Spectrum estimation for a portion of sphere. 



Hi. 1 Estimator of the Ci spectrum. 

The temperature field is measured on a piece A of the sphere of solid angle tt obs . We 
will study the following estimator for the Cf 

1 J 

a 



^yEi^i 2 ( IIL1 °) 

m=— I 

where the a; m are computed on tt obs : 

ai m = J T(tt)xY™(tt)dtt (IIL11) 



The Ci will therefore represent an estimation of the spectrum of a temperature map which 
is set to zero outside the observed region. We define the function W(tt) on the sphere 
such that W = 1 on A and anywhere else. 
The temperature field under study is: T'(tt) = (T x W)(tt). 

Hi. 2 Computation of the bias of the estimator. 

The function W can be expanded on spherical harmonics basis as: 

oo I 



= EE h imYrm (ni.12) 

1=0 m=-l 

hm = J W{tt) Y^(tt) dtt = J Y™(tt) dtt (111.13) 

S 2 &obs 

Then the a; m are: 

aim = I T{tt)W{tt) Y^itt) dtt = E a l'm'k"m" [ Yp'Yfi'Y™ dtt 

J II / 111 // J 



S 2 l',m',l",m" s2 



Using the property Y™ = (-l) m Yf m : 



aim = Yl a l'rn'k»m»{-l) m f Y™' Y~ m Y{f dtt 
V,m',l",m" £ 2 



The integral over the sphere of the product of three spherical harmonics is ( |Brink et al, 1962] 
and |Messiah, 1964] ): 



« - f° +i n: 1)(2c+1) (: \ ;) (; i s) 



s2 
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where ( „ 1 are the Wigner 3j symbols related to the Clebsch-Gordan coefficients 



by: 



a (3 7 



<aba(3\c- 1 >= + \ f) 



Recall that these coefficients are real. 



a>im — Yli tt i'm'W"(~l) 



m / (2l+l)(2l'+l)(2l"+l) 



47T 

I'm'M'.m" 



0-lm — ^2 0-L'M'i>L"M"{ — ^Y 



m /(2/+l)(2L'+l)(2L"+l) 



Ait 

L'M'L",M" 



'"„) 




I 


I" 


m") 


Co' 








I 








—m 


M" 


)( 


o 



*2/m I O'lmO'lm 



2/ + 1 

5^ (a,/ m /aZ^T) (b v , m „bv>M») x V(2/' + 1)(2Z" + 1)(2L' + 1)(2L" + 1] 

m! -m m")\M' -m M" ) [ ) { 



x 



where the sum runs on l',m , I", m , L', M , L , M . 

The aim are random variables. We want to compute the average effect of finite size 
maps, therefore we shall perform the average over sky realisations (ensemble averag^l). 
Recalling equation ( 111.31) . one obtains: 

1 1 

<Cl> = -^Yl J2 C v b v , m „b L „ M ,,(2l' + 1) y/(2l» + 1)(2L" + 1) 

m=-l l',m',l",m",L",M" 

(V I l"\(V I L" \ (V I l"\ (V I L" 
X \rri -m m") \m' -m M" ) ^0 0J \0 

Now, remark that the coefficients of the Wigner 3j symbol product do not depend on m 
nor m', therefore one can use the orthogonality relations 



a/3 

Defining: 
we obtain: 



^ v / (2c+l)(2c' + l)[ 



a b c\ ( a b c'\ . . 
a (3 7 J U P i) = cc n ' 



i 

i 



J— J ^ b l"m"bl"m" (III. 14) 



21 

m"=-l" 



<Ci> = j-J2Ci'B vt (2l'+l)(2l" + l)(l g qY (IIL15) 



2 In the following, we will drop the subscript for ensemble average notations: < ■ • ■ > ens < • • • > 
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The expression of the 3j symbols can be found for instance in Brink et al, 19 62 . We 
found recently that this calculation has already been published in Hivon et al, 2 002 . 
The above equation provides the relation between the true C\ and the average < Ci > 
over many sky realisations of the G% coefficients measured over a patch A of any shape of 
the sphere. This can be written in matrix form: 



<a> 



[M) lv x [C,,] with M IV = ^ + i^ (2r + J qY(HI.16) 

l"=0 ^ ' 



This result is valid for any weighting function W on the sphere. The relation between the 
C\ and the C\ is linear. The matrix Mw depends only on the map shape and not on the 
initial Cy spectrum. Permutations of 2 columns of the 3j symbols change their values by 
a phase QMessiah, 1964] ) so (2/ + l)M w = (2V + l)M n . The M w elements are positive 
or null. 

When the full sky is observed, By = B 6y and the 3j symbols imply I = V whence 
Mw is the identity matrix. 

On the other hand, when observing a portion of the sky, each measured C\ is a mixture 
of the Ci whose weights are given by the matrix Mw- 

The properties of the matrix are illustrated in figure HI The bottom graphs show 
line / = 1500 of matrix Mw for a 10 x 10 deg 2 map and lines I = 500, 2400 where the 
peaks have been shifted to I' = 1500 for easy comparison: the width of the peak is nearly 
constant with I, that is to say that the matrix is nearly band-diagonal. 
The top graphs show the Mw element values of line I = 1500 for maps of various size 
and shapes. As the map size decreases, the peak width gets larger and a given Ci mea- 
surement involves a larger band in Cy. The width of the peak is inversely proportionnal 
to the caracteristic width of the map and it does not depend very much on the shape 
provided it is not excessively elongated. This is somewhat equivalent to what we have in 
Fourier analysis on small intervals. 

The C\ are coefficients in the correlation function expansion ( 111.51) . The P\ are orthogonal 
on [0, 7r] and form a basis on this interval. On a partial map of characteristic size 9 C these 
functions are "nearly orthogonal" if they differ by more than 51 ~ tt/6 c : loosely speaking 
if they have not the same number of roots in that interval. Therefore, for a given corre- 
lation function, the C\ will be mixtures of the "true" C\ over a range ~ 51. 

If Q b s is the solid angle corresponding to the observed patch of the sky, using the 
completude relation for the spherical harmonics, we have for the Bi normalisation: 



oo oo I „ 

^2(2i +i)Bi = Y. bi ™^ = / \ w $) i 2 d ' n 

1=0 1=0 m=-l 



S 2 

Until the end of this section, we assume that W = 1 on A and zero elsewhere. Thus we 
have: 

oo 

^2(2l + l)Bi = Q obs (IIL17) 

1=0 

The C/ are biased relative to the C\ since one observes only a fraction of the sphere, the 
rest being 0. Assuming that the function T is isotropic, to compare the C\ with C\ one 
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has to take into account the sky coverage normalisation: 

~ 47T ~ A.7T 

< Cf > = — - x < Ci > = M$ x C v where Af$ = -— M IV (III. 18) 

In £/ie following, when talking about comparison between C\ and C\, we will assume that 
the Ci have been renormalised as describe above. 

Using the 3j orthogonality relation ( |Messiah, 1964| , appendix C relation 156) which 

gives $2(21' + l)(nnn) =1 and ec l uation pTT7j) we get J2 M w = 1 - 

l' V u „ / v 

Therefore for a given < Of* >, the sum of the weights of the Gy is equal to one. 



The Bi can be obtain using usual computer technics (see Gorski et al, 2005] ) . Then 



it is easy and fast to compute numerically the matrix Mw and check the systematic effect 
of finite size map. 

The figure [2] shows the B\ spectrum and the M\y matrice for a 10 x 10 deg 2 square map. 

Hi. 3 Computation of the covariance of the C\ 
estimator for spherical cap maps. 

We have tried to compute the variance of the Ci estimator following the same line 
of calculation. Although the result does not look like too complicated, it is of little use, 
because it would cost a lot of computer time to be numerically calculated in the general 
case. However it can be calculated in the simple case of a spherical cap domain. 

If, in equation ( 1III.11I) . we replace the temperature field by its expression (111.11) . we 

get: 

l\m\ 

where we have set: 



B hmilm = J w(n)Y^(n)Yr(3) dQ 



Averaging over sky realisations and using equation (III. 31) . one gets (see details in ap- 
pendix [B]) : 

< Ci > = ^ 1 - ^2 c v I B imi'm' I 2 for -I < m < +1 , < I' < oo , -f < m' < +1' 

ml'm' 

and provided that the ai m are gaussian random fields: 

_ _ 2 +1 +L +v 

< CiC L > - < Ci >< C L > = — 2J I ^ B Vm ,i m B Vm i LM | 2 

^ ' m =-l M=-L V m'=-V 

The last two results are valid for any weighting function W on the sphere. W could 
include an additionnal weighting, for instance to treat border effects. The results of the 
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former section can also be adapted to that case. Until the end of this section, we will 
assume that W = 1 on A and zero elsewhere. 

This equation simplifies a lot if one considers a spherical cap centered at the north 
pole with border at polar angle 9 C . Using Y l m (9,4>) = A™(cos(0)) e 1 " 1 ^ (where Xf 1 are 
the normalised associated Legendre polynomials), we get in this case: 



Bi mL M = 2tt 5 mM B lLm wherd B lLm = / Ap(cos(0)) A£(cos(0)) dcos(0) 



Replacing in the general expression for the covariance, one obtains: 



< a > = |^ J2 E ( fi ^) 2 ( IIL19 ) 

m=—l l'>\m\ 



and for the covariance: 



< QCl >-<d><C L >= {2l + ypL + l) E f E °v Bii'mBu.m ] (HI.20) 

The Bwm coefficients could be computed by numerical integration but this would be very 
slow at high /. The calculations were done using recurrence relations. 
The figure [3] shows the correlation matrix computed from the covariance formula above 
for a 10 x 10 deg 2 spherical cap. 

Recurrence relations 

In order to derive the following recurrence relations for the Bu/ m coefficients, we have 
used the usual recurrence relations among the normalised associated Legendre functions 

n c 

( |Gradshtein et al, 1980| ) which imply integrals of the forrro J xX™(x)X™(x)dx. The 

i 

c 

latter integration can be in turn related to integrals of the type J A™(x)A™(x)<ir by 

i 

c 

evaluating j(l — a; 2 ) (A™' A™" + X^'X^^dx and using the Legendre differential equation 
i 

(where ' and " means first and second order derivations with respect to x). One gets: 



Bii+i )m — — cA 



r(c)x T+1 (c) + y^gp^j (^|±3 mc)? + ^ 3[ x r+1 (c)f 



3 We have B Wm = B Vim = B W - m G K , B oo = 37 (1 - cos(0 c )) , B Wm = for | m |> min(/, I'). 
From the relation E>T(^)^ m (^) = ^Pi^h^h) we get £ B Um = ^ n obs 

m m 

4 We set c = cos(# c ) and x = cos(#) 



12 



And for | k — I |> 2, setting b = (7 + l)(k — l){k — m), one has: 
b x B ikm 



k 



2k + 1 21 + 1 k-m l + l-m 



2k-l\l2l + 3\lk + m\ll + l + m 
-(k-m)l(l + l)cAP(c) - 



(/c-/-2)(/ + l + m) S,+i,fc-i, 



2/ + 1 l + l+m 



21 + 3\ l + l-m 



2k + 1 / k-m 

2k^i\lkT^ Xk - l{ c) 



(i + l-«i(c) AHc) 



(III.2i; 



x I fc y|^v l + l + ^ ( / + 1 + ™) cA m(c) - (*(* + 1) " m 2 )XT(c) 
This heavy recurrence relation has been verified to be numerically stable up to I ~ 2500. 



The only numerical integral to perform is Bu m . As A[ 



[x 



[1-x 1 



Q{l — m) (where 



Q(n) stands for a polynome of degree n), the integrand of Bu m is a polynome of degree 21 
and it may be exactly computed using the Gauss-Legendre integration method with I + 1 
points. All other B^m can be computed using the recurrence relations flHI21j) . 



With similar technics Wandelt et al, 2001 derived different equivalent recurrence rela- 
tions. 



Approximation for small maps 

For small 9 C (8 C < 0.15), one can make the approximation ( |Gradshtein et al, 19 80 1 ) : 



Xl (C ° SW) ^ ~F~y 47r(Z-m)! Jm{W) 
where the J m are the cylindrical Bessel functions. Therefore: 



V 4vr 4tt (I - m)\ (h - m)\ lfl m J u; 1 ; 
v 

where the last integral is a Lommel integral ( |Gradshtein et al, 19 80 1 ) . Then: 



B, 



l\lm — 



'21 + 1 2h + 1 (l + m)\ {h +m)\ 1 9 C 



Att 4tt (l-m)\ (Zi-m)! l?l m If- 1 2 

x {u m {i x e c )j' m {ie c ) - hJ m (w c )j' m (he c )} (111.22) 

where J' m (z) stands for the derivative of J m {z) relative to z. 

For numerical purposes it is easier to replace {U m (li8 c )J' m (W c ) — liJ m {l0 c ) J' m {li6 c )} by 
{hJ m (W c )J m+ i(li9 c ) — lJm(li9 c )J m+ i(W c )}- Figure M compares the variance of the C/ 2 
computed with relation (1III.21I) with the calculation done with the above approximation: 
the agreement is excellent for f sky = < 6 10~ 3 (9 C < 0.15 rd). 



Behaviour of the variance as a function of 9 C 

If 19 c < m, Jm(W c ) behaves as (W c ) m which gives a negligible contribution if m is large. On 
the other hand, if m < 19 c one can use the following approximation ( |Gradshtein et al, 19 80 1 ) : 

[2 1 

J m (W c ) ~ \ = cos(/6' c + m ) 

V 7T y/W c 
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where m is a phase depending on m and slowly varying with 19 C : <\> m — > — m| — | for 
W c — > oo. 

In that approximation, from equation (1III.22P one obtains^): 

B^im ~ — f 1 — - 6' c = 6 C + corrections depending on m (III. 23) 

One can use this approximation in equation (1III.20P in which we set I = L for the variance 
calculation. 

If only the central peak of Bf lim as a function of l\ is considered — I |< n/9 c ) the 

partial sum ^2 in equation ( 1III20I) behaves as 9 C because Bf lim is a function with a peak 
V 

height proportionnal to 9 C and width ir/9 c . A more systematic calculation shows that, 
due to the high values of Ci at low I, one can not neglect the contribution of Bf lim at the 
left of the peak (li < I — m/9 c ). 

The sum over m in equation 20[) introduces a factor ~ 219 C . 

If one considers only the central peak (and using the normalisation in formula (1III.18P ). 
one gets a 2 ~ R ~ — - , so <jq r would scale as f sk y. 
Numerically we found that a good approximation for small 9 C is: 

g c ; ( 47r ) mioA\ 
a c* ~ rr- nr ( IIL24 ) 

V JskysJ q c 

For large I values, the peak wing corrections lead to corrections of higher order in 1/9 C : 

a 2 r (4m) b(l) 



Hi. 4 Discussion and comparison with simulations 

In the following the C\ coefficients are computed from simulated CMB maps using the 
well known Fast Fourier Transform (FFT) technics ( |Natoli et al, 1997| ). 

Simulations have been realised using the HEALPIX package ( |G6rski et al, 20 05 1 ) . 
Pure CMB maps of various sizes and shapes have been simulated using the program 
SYNFAST with l max = 2500 for the standard ACDM spectrum. The temperature field 
outside the various observed regions is set to zero. 

The pixel size was chosen small enough (NSIDE = 2048) in order not to bias the 
Ci reconstruction up to l max = 2500. The pixel smoothing windows has been chosen to 
be that of a NSIDE = 8192. So the value in a pixel is very nearly the value of the 
temperature field at the center of the pixel. No detector noise nor telescope lobe effect 
was added. 

The Ci were reconstructed using the program ANAFAST up to l max = 2500. That is 
to say that ANAFAST computes our C\ as defined in equation flIII.lOjl . 

Using full sphere simulations we checked that the reconstructed C\ were in perfect 
agreement, even at large /, with the generated ones. Table (CQ) shows the generations used 
in the following discussion. In the following f s k y = is the fraction of the observed sky. 

5 This equation and (|III.19p explain qualitatively the behaviour of the matrix Mw . 
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map size 


f sky 


number of generations 


47r sr 


1 


579 


20 x 32 


1.54 10" 2 


228 


17 x 17 


0.697 10- 2 


579 


R = 9.585 


- 


2279 


10 x 10 


0.242 10~ 2 


228 


6x6 


0.872 10- 3 


228 


4x9 




228 


i2 = 3.385 




228 



Table 1: Description of simulations: R = 6 C means a spherical cap with polar angle 
9 C deg, other maps are labeled A9 x A<f> deg 2 . 



Figure shows that the bias computed using formulae (IIII.16|) and flIII.18|) is in excel- 



lent agreement with the simulations. Figure [U] illustrates the fact, with 17 x lldeg 2 maps, 
that this systematic bias may not be small with respect to the variance. In the present 
figure it is of the order of 2a for I > 2000. We note also that the correction can be im- 
portant at lower I although dominated by the variance. The Cf 1 values are systematically 
above the input C\ when 200 < / and are oscillating with a period 81 ~ tt/6 c . The matrix 
Mui (1HI.16P has long oscillating quasi-symetric tails (see the figure HJ). The input C\ spec- 
trum used is a rapidly decreasing function. Therefore the left tail contribution is much 
larger than the right tail one, making the bias always positive at high I. The observed 
oscillations are a consequence of the Mu> matrix oscillating tails and of the global shape 
of the Ci spectrum, but are not, at first order, the image of the acoustic peaks. A smooth 
input Ci spectrum without peaks leads also to an oscillating CP spectrum with only a 
sligthly different pattern. For a more rapidly decreasing input spectrum, the amplitude 
of the oscillations increases. 

The variance has been computed numerically using ( 1III.20j) and the recurrence rela- 



tions flIII.211) . The results are compared with various simulations in figures [7] and [HI They 
show that the variance can be predicted accurately independently of the map shape, al- 
though the calculations were performed for spherical caps. The figure \E\ shows in detail 
the agreement between spherical cap map simulations and calculations. We can even go 
further and compare the correlation matrix obtained from simulations to the calculations 
for spherical caps: 

„ fl rN < Ci C L > - < Ci >< C L > 

Cor(l,L) = 



<C 2 >-<Q > 2 v /< C\ > - < C L > 2 

We will do the comparison for a few lines of the correlation matrix. Figure [TU] compares 
computation and simulation for two lines (I = 2001 and I = 1500) of the correlation 
matrix around the correlation peak for f s k y = 0.697 10~ 2 (see table (CQ)). Top plots 
are for a spherical cap (2279 simulations) and bottom ones are for a square map (579 
simulations). When I, L are both larges, the dispersion is rather small. The top left plot 
shows the very good agreement between computation and simulations. The bottom left 
plot shows that, even for a square map, the correlation peak is well reproduced. Outside 
the peak, the correlation level is fairly good but out of phase due to the map shape. The 
same conclusions can be drawn for the right plots although dispersion in the simulations 
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is larger because of the large variance we have for small values of /, L. In conclusion, the 
correlation peak shape is rather independent of the map shape. We note also that the 
correlation matrix outside the diagonal region has a periodic structure of period ~ 2tt/9 c 
which is a consequence of formula (1III.23I) . 

Figure [H] shows crp} R /aci (47r) versus f s k y for various I. The variance (cx~ fl ) does not scale 

as f~ k y as often assumed (see the discussion in ( |Tegmark, 1997| )). The black dashed 
curve shows cr^ R / 'crc,(47r) = i — . 

For / ~ 500 the variance scale as f sk y as expected from the lowest order calculation 
(formula (1111.241) ) represented by the solid black curve. For larger I values, we have 
to take into account the correction terms of formula flIII.25p . This indicates that the 



dependance of the variance on the map size depends also on the Ci spectrum shape. 
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Figure 2: The top picture shows the Bi spectrum for a 10 x 10 deg 2 square map. 
The bottom picture shows the matrix Mu> for the same map (The LUT as been optimised 
for readability purpose). Axes range from upper corner left (/ = 0,1' — 0) to bottom right 
corner (Z = 2500, V = 2500). 
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Figure 3: Correlation matrix for the C\ a 10 x 10 deg 2 spherical cap. 

The top left corner is for (I — 0,1' — 0) and the bottom left corner is for (/ = 2500, 

2500). The LUT ranges from -1 (dark) to +1 (white). 
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Figure 4: The top picture shows line / = 1500 of matrix Mu> for I' ranging from 1100 to 
1900 and for maps of various sizes and shapes. The bottom one shows line I = 1500 of 
matrix Mu> for a 10 x 10 deg 2 map and lines I = 500, 2400 where the peaks have been 
shifted to /' = 1500 for easy comparison. 
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Figure 5: Systematic effect on C\ for various size maps. The thick black curve represents 
the Ci coefficients used in the simulations. The colored curves show the average < Cj 1 > 
reconstructed from simulations for various map sizes. The thin black curves superimposed 
on colored ones show the results predicted by formula (lIII.16p . 
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Figure 6: Finite size map systematic effect on C\ for a 17 x 17 deg 2 map. The thick 
black curve represents the C\ used in the simulations. The red curve show the average 
< CP > reconstructed from simulations. The two pink curves show the one sigma de- 
viation reconstructed from simulations. The blue curve shows the results predicted by 
formula flHI16j) . 
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multipole moment 

Figure 7: Finite size map effect on C\ variance for various maps. The thick black curve 
represents the variance of the G\ for the full sphere. The colored curves show the variance 
of the CP reconstructed from simulations. The thin black curves superimposed on colored 
ones show the results predicted by formula 21 [) . 
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Figure 8: Finite size map effect on C\ variance for various ~ 36 deg 2 maps with dif- 
ferent shapes. The zoom compares the result of simulations for spherical cap with the 
formula <jlll.21j> . 
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fsky = 1 

flower) fsky=1 .543 1 0" 2 with Leg end re Pol. 
(lower) with Bess el tunc. 

fsky=6.979 10" 3 with Legendre Pol. 

with Bessel tunc. 
fsky=2.421 10" 3 with Legendre Pol. 

with Bessel tunc, 
(upper) fsky=a.723 10" 4 with Legendre Pol. 
(upper) with Bessel tunc. 
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Figure 9: Validity of the small polar cap map approximation. Comparison of the vari- 
ance computed with Legendre polynomials and Bessel function approximations (see sec- 
tion QUOD ). 
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Figure 10: Comparison of the computed correlation matrix with simulations for 
17 x 17 deg 2 maps. Black curves correspond to the computed correlation for a spher- 
ical cap. Red curves represent the correlation computed from simulations: the top figures 
correspond to a spherical cap and the bottom figures to a square map. The lines I = 2001 
and I = 1500 of the correlation matrix Cor(l, L) are drawn versus columns L around the 
correlation peak. 
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Figure 11: Square root of the variance of C\ relative to the full sphere variance versus 

_ i 

map aperture 6 C for various I. The black dashed curve represents the simple f sk l scaling. 
The black solid curve represents the first order scaling (equation (1111.240 ). 



26 



IV Angular correlation function estimation for 
a portion of sphere. 

iv. 1 Estimator of the angular correlation function. 

The angular correlation function £(7) is estimated as in section lii.3l but with the 
integration limited to the portion of sphere A. Let's define : 

W A (Vt) = 1 if Q e A and elsewhere, or W(Vt) > in A and elsewhere. 

The angular correlation function can be computed without ponderation on the tem- 
perature field. We define: 

ffr) = Jfh) / ^1^2^)^2)5(01.^-008(7)) 



Ax A 



AM (7) 

S 2 xS 2 



/ dn 1( in 2 T(n 1 )w A (n 1 )T(n 2 )w A (n 2 ) 5(Oi.o 2 - cosfr)) 



or more generally with a ponderation W(Q): 

f(T) = I dn 1 dn 2 T(n 1 )W(^T(Q 2 )W(Q 2 )5(Q 1 n2-cos^)) 



AxA 



= j^t-, f dnidn 2 r(n 1 )w(ni)T(n 2 )w(n 2 ) < y(n 1 .n 2 -cos(7)) 

S 2 xS 2 

with 

N A ( 1 ) = J dfiitZfia - cos(7)) 



/ df2idf2 2 W^(fti) W^(0 2 ) 5(fii.fi 2 - cos(7)) 

5 2 x5 2 



We have T(Q) = aim ^(Q) and we define: 



LAI 



W(Q) = £& LM F L M (0) 

LM 

w A (n) = Eb A M y L M m 



LM 



The above integral becomes: 
A/" A (7) £(7) = £ W^b LlMl a hm2 b L2M2 



x / dflidfla 5(^.^2 - cos(t)) F^^i)^ 1 ^!) ^ 

S 2 xS 2 

where the sum goes over {li, mi, Li, Mi, / 2 , m 2 , L 2 , M 2 }. 
The computation is described in details in appendix O We get: 

2i + 1 lEet + DBLftM 



<«?>>= E^QflW 



4' ' / E (2i + 1) Si 4 Pl(7) 

L 
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where Bi and Bf are defined respectively with bi m and b A m according to (lIII.14p . and 

A/^(7) = 2vr^(2L + l)Si 4) P L (7) 



For 7 > Q m ax, where Q ma x is the largest angular distance on the map A, £(7) is not 
computable. 

• If W = W A , the temperature field is not ponderated on A and the estimator is unbiased: 

~ 2/ -I- 1 

<£(7)>= E-^^W = 

• If W ^ W A , the temperature field is ponderated on A by the positive function W(Q) 
and the angular correlation function is biased. 

• If the temperature field is ponderated on A by the positive function and we 
replace J\f A ('j) by: 



(7) = J dn t dn 2 w(fii) w(n 2 ) s^.^Z - cos( 7 )) = 2n ( 2L + b l p M 



S 2 xS 2 



the estimator is unbiased. 

Note that J\f A ('j) is proportional to the number of pairs of directions separated by an 
angle 7 that can be done on A. 

iv. 2 Comparison of the estimator with the < C\ > 

< £(7) > can also be expressed relative to the ensemble average < Ci > of the bias Ci 
computed in section lin~2l 
We have (c = cos(7)): 

N A (c)£(c) = J2 a^ai>m> J dQdtt' 5{Q .Q' - c) Yp(fi) Yf {M) 

Iml'm' AxA 

The integral is computed in appendix [Dl Computing the ensemble average gives: 

mm' X (^IttPl^C) B^Mlm Bj^Ml'm'^) 

Iml'm'LM 

= 27T £ QP L ( C ) I B LMbn | 2 
ImLM 

= 2tt£ P l (c) £ Q I B iJWm | 2 
= 27r£P L (c)(2L + l)<C L > 



We finally obtain: 

Af A ( 7 ) 
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<£( 7 )>= 2^(2/ + l)^g) <a> 



iv.3 Covariance of the estimator of the angular 
correlation function on partial map. 

Let's define C\ = cos (71) and c 2 = cos (72). 
M A { Cl )l{ Cl ) = E^Oiimi / dn 1 dQ! 1 5(ShS l -c 1 )Y^{n 1 )Y^' l {Q! 1 ) 

AxA 1 

where the sum S oes over 

M A { Cl )M A {c 2 ) < a Cl )l{c 2 ) > 

= < G^mi fliimi a l' 2 m' 2 a l 2 m 2 > J dVtidVt^ 5(^1.^ — Ci) 

x / dQ 2 dQ' 2 5(tt 2 .tt 2 - c 2 ) Y^fii) Y^'^'i) F™' 2 (fi' 2 ) 

AxA 1 2 

If the temperature field is real and gaussian, we can compute the fourth moment of the 
ai m . That leads to 3 terms: 



1- / c v c V2 ■ -y,r i (fl > 1 ) ^(fl 1 ) ^(n' 2 ) ^(n 2 ) 

2- / a i a---^k^)^(fii)^(^)^Hfi 2 ) 

3 _ / (-l^i-Hmcijc?,, • • (ni) V 1 (no V*' 1 to V 1 M 

The first term gives < A/" A (ci)f(ci) >< N A (c 2 )l(c 2 ) >. Using Y^ m = (-l) m Yf m and 
playing with indices, it is easy to demonstrate that the second and third terms are equal. 
We have: 

N A { Cl )N A {c 2 ) (< l{ Cl )l{c 2 ) >- <£( Cl ) >< £(c 2 ) >) 

= ^( Cl )^(c 2 )y(e( Cl ),e(c2)) 

= 2 x £ 

x / dQidfii <J(fiI.fi{ - ci) Y^i^Y^i^) 

AxA 1 

x / dn 2 dn' 2 8{n 2 M 2 - C2 )Y™'{n 2 )Y^{n' 2 ) 

AxA 1 

= 2x £ ^W^i) & 2 ) 

where: 

i\r m ,{c) = J dmn'5{nxi'-c)Yf{a)Yr{n) = 7£^) 

AxA 
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is computed in appendix [D] Going back to the computation of the £(c) covariance: 

N A { Cl )U A {c 2 ) V{l{ Cl )^{c 2 )) 

= 2(2vr) 2 E ^Qj^^fite) 

= 2(2n) 2 E E C^C^ Pl(ci)Pl'(c 2 ) Blmi^Blmi^ B^Mn^BLtM'hnu 

hmil^m^ LML'M' 

= 2(2tt) 2 E ^(ci)Pl'(c 2 ) E B LMhmi B L , M , hmi E ^ ^ 

LML'M' \hmi J \V x m' x 

= 2(2tt) 2 E Pl{ci)P L '{c 2 ) I E B L/M 'i m B L ~ M ^ | 2 

LML'M' Ira 

= 2(2tt) 2 E Pl(ci)P^(c2) ( E EQW^I 2 ) 

The term in parenthesis can be expressed relative to the covariance Vq(L,L') of the C; 
(see appendix [Bj: 

A^( Cl )A^(c 2 ) l/(e(ci),e(c 2 )) = 2(2tt) 2 E Pl( Ci )P l ,(c 2 ) ( ( 2L+1 )f +1 ) y g (L, L')) 

LV v 7 

= (2tt) 2 E (2L+ 1)(2L' + 1) V 5 (L,L') P L {^)Pv{c 2 ) 

LV 

Using the value of Af A ( r y) computed in section EH one finally obtain for a partial map of 
any shape: 



n£(ci),e(c 2 )) = < e(ci)£(c 2 ) > - < e(ci) >< e( C2 ) > 

( 27 r) 2 £ (2i+l)(2i'+l) Vg «/) PK7i)Pj/(7a) 

— Ar^( 7l )Ar^( 72 ) 

E (a+i)(2i'+i) v~ ^(71)^(72) 
w _ 

(E (2i+l) P L ( 7 i)) (E (2i'+l) B£ P L /(72)) 
L r,' 



Recall that the Vg ; (7, /') can only be computed easily if the partial map is a spherical 

cap. The figure [T2l shows the correlation matrix of the £(7) for a spherical cap of area 
17 x 17 deg 2 : one sees that the level of correlation is very high. 
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Figure 12: Correlation matrix of the £(7) for a 17 x 17 deg 2 spherical cap. The top left 
corner is for (7 = 0,7' = 0) and the bottom left corner is for (7 = 19.1,7' = 19.1 deg). 
The LUT ranges from 0.5 (dark) to +1 (white). 
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V C\ estimation using the angular correlation function 



for a portion of sphere. 

v. 1 Using integration of the angular correlation 
function. 

Suppose that we measure the angular correlation function £(7). 
For the full sphere we have: 



Q = 2tt J £(7) P(cos( 7 )) dcos( 7 ) 
7=0 

On a portion of sphere, we may have an estimation of the C\ by performing the integration 



Cf = 2vr J £(7) P z (cos(7)) dcos( 7 ) (V.26) 
7=0 

where 9u m is the maximum separation angle obtainable for that map. 
We demonstrated that £(7) is unbiased, so: 



00 

< £(7) > = £(7) = ^ £ (21 + 1) Q P(cos( 7 )) 

Thus we compute: 



Aix 
1=0 



< Cf > = 2n J <£(7)> Pi(cos(7)) <icos(7) 
7=0 

@lim 

= 2ir J £(7) P(cos(7)) dcos(7) 
7=0 

OO Olim 

= 2vr^ E {2l' + l)C v j P r (cos( 7 )) P,(cos( T )) dcos( T ) 

l'=Q 7=0 

Using the definition of the Bu m , the Pi and the A° we have: 



< Cf > = 2n £ Jm C V B, 



21+T ^i'lo 

i'=0 v 

with P// ro = Bifi (du m ). 

H As we have a sharp cut-off at 9u m , the obtained spectrum oscillates strongly around 
the theorical value (that is somewhat equivalent to the Gibbs phenomena we have with 
the Fourier transform), and that method is not usable. The figure [TH1 shows the predicted 
Cf for a spherical cap of angular aperture 9u m = 19 deg. The level of oscillations is very 
high. The "wave length" of the oscillations is about Al ~ This problem is discussed 
in the next section. 
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One can compute the covariance of the Cf: 

<Cfdf,>= (2tt) 2 J <l{ 1 )l{ 1 ')> P ^ (cos(7))P^(cos(7 , ))rfcos(7)rfcos(7 , ) 



Using the variance of the £ computed in section liv.3[ we obtain: 

V(Cf,C*) = J2 M ih V dl {kA) M l% (V.27) 



with 

M lh = (2vr) 2 (2h + 1) J ^1 fl( 7 ) fll ( 7 )dcos(7) (V.28) 

o 

t;.^ Using integration of a smooth apodization of 
the angular correlation function. 

In order to avoid the problem shown in figure [13] due to the sharp cut-off at 9u m , the 
correlation function £(7) is multiplied by a function going smoothly to zero at 6u m . 

A usefull function is: F tan h(0) = ~(1 — tanh(^ a )) where 6q is the cutting angle 
and A caracterises the width of the smoothing. 
One could also take F erf (6) = ±(1 - erf(^)) . 
The estimator in (IV. 261) is changed to: 

Cf = 2?r y £( 7 ) F( 7 ) Pi( cos (7)) dcos(j) (V.29) 
7=0 

In that case the average value of the estimator is: 

<Cf> = 2tt / < £(7) > F( 7 ) P ; (cos(7)) rfcos(7) 
7=0 

OO Slim 

= §£(2J' + 1)CV J ^(cos( 7 )) ^(cos( 7 )) F( 7 ) dcos( 7 ) 

i'=0 7=0 

It has to be computed numerically. 

The covariance keeps the same form as (IV. 271) if we replace Mn t in the former section by: 

M lh = (2vr) 2 (2h + 1) J ^l^p( 7 )P ;i ( 7 )F(7)rfcos(7) (V.30) 



9¥ Apodization effect for small maps. 
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We compute: I = f P//(cos(#)) Pz(cos(6*)) F{9) dcos(9) . For 8u m small, we may 

6»=0 

approximate Pi(cos(9)) ~ Jo(W) and we have (sin(6>) ~ 9): 



I = J J {1'0) J (W) F(9) 9 d9 



o 

e 

We define: L(9) = J J (l'9') Jo(W) 9' d9' which is a Lommel integral. Integrating / 

o 

by parts: 

I = [L(9) F(9)]l Hm ~ J L{6') ^F{0') dff 

o 

But L(0) = and we set the smoothing function to be null at 9u m i.e. F(6u m ) = 0. Thus 
[. . . ] i™ = o and we have: 



To perform analytical computations, we take: F{9) = F er j(9) = |(1 — er/(^ a )) , we 
have: ^§F er f = — exp (— 1 and we obtain: 



Jl{6) ex P (-^^) d9 
o 



The Lommel integral L(9) is ( |Gradshtein et al, 19 80|): 

L ^ = i^rp WWm - i'MW)J'S'o)\ 

where Jq(x) is the derivative with respect to x, and we make the approximation (for 
W Um > 1): J (W) ~ cos(Z0 + <p) with <p -> -(n + |)f if 9 -> oo . For J , this 

remains a good approximation even for values as low as 19 ~ 4. 
We obtain: L(9) = [^M^M + sin((r+;) e+2 ^ 

The second term is of order with respect to the first one and will be neglected. We 
finally have: 

i <* *foT ^ (-^) * 



o 



JW J I'-l eXP I A2 J Qt/ 



where we have replaced the integration limits because A ^ 9q so that the integrand is 
null outside the limits. Setting z = 9 — 9 we have: 

sin((f - 1)9) = sin((Z' - /) (0„ + z)) = sin((/' - Z)0 O ) cos((Z' - l)z) + cos((/' - Z)0„) sin((Z' - l)z) 
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Thus: 

f +00 

1 « a7!F "^fer^ / «»((?- Q*)exp(-£)d* 
L —00 

+00 

+ cos((J / -Q0 o ) / ^fe^exp(-fi)^ 

—00 

The second integral is null by symetry: 



rOG 



I oc - / cos((r-/),)exp(-^)^ 



—00 

00 

As / cos((Z' - l)x) exp(-a; 2 /A 2 ) dx = y/n A exp(— |(Z' - Z) 2 A 2 ) , we have: 



#0 sin((Z' - Z)# ) , 1,,, /\2 a 2\ 



and finally: 



(2f + 1) c r A= 5^=^ exp(-i ( r - () 2 A 2 ) 

The value of a < Cf > is given by a kind of convolution: 

• The width of the sine goes like ^, so as #0 goes to zero, the width will become large. 

• The width of the gaussian goes like 4- (cr = so the larger the width of the 
smoothing function (large A), the smaller the width of the gaussian in the last 
formula. 

The oscillations described in the previous section (see figure Q3]) are due to the contribu- 
tion of the low multipoles at high I: if there is no apodization, as there is a lot of power 
in low I multipoles, the high I multipoles will get a lot of power transfered from the low 
ones and the oscillations will be large: 

- If the apodization is too sharp, the width of the gaussian will be very large and the 
"convolution" will be dominated by the sine. The "sine" decrease slowly (1//), the trans- 
fert of power from the low multipole to the high ones will be important and so will be the 
oscillations. We need a smooth apodization function. 

- If 6q is large the width of the sine will be small and the oscillation of the sine will remain 
imprinted on the spectrum. We must take 9 not near 6u m . 



The figure dH shows the effect of various apodization parameters on the predicted C\ 
(formula IV. 291) . The apodizations are applied to the "theoretical angular correlation 
function" with 9u m ~ 19 deg corresponding to 17 x 17 deg 2 spherical cap maps. 

• for #0 — 15 deg, A = 0.5 deg: the spectrum oscillates at low I due to the sharpness 
of the filter. At high I it is undistinguishable from the theorical one. 



35 



• for 6q = 10 deg, A = 5 deg: the angular correlation spectrum is not apodized enough 
near 9u m and, as explained above, the low multipoles C\ values correlate with the 
high multipole ones leading to oscillations at high I in the reconstructed spectrum. 

• $o — l deg, A = 2 deg shows that one cannot lower too much 9 without loosing 
information. The reconstructed spectrum becomes highly biaised and smoothed. 

• #o = 10 deg, A = 2 deg shows an example of good apodization. The reconstructed 
spectrum differs only at very low multipole values where the size of the map becomes 
too small. 

Using IV.27I with IV. 30} the top plot of figure [15] shows the variance 1(1 + l)ag^/2ii for 
(9 = 10 deg, A = 2 deg) and (9 = 10 deg, A = 5 deg): because of the poor apodization 
of the second filter, the variance increases a lot at high multipoles. 

The bottom plot of figure [T5] shows the correlation function for I = 1500 for (9 = 
10 deg, A = 2 deg) and (9 = 5 deg, A = 2 deg): the width of the correlation peak 
depends on the value of 6q and is greater for the second filter as the accessible separation 
angle in the angular correlation spectrum is lower in that case. 

v. 3 Practical C\ reconstruction from the angular 
correlation spectrum. 

The above analytical results were applied to a set of 173 simulated 17 x 17 deg 2 square 
maps to check the quality of the C\ reconstruction. Computing the correlation spectrum 
is time consuming, therefore, for practical purposes, the temperature field has been cal- 
culated on a rectangular grid with constant 9 and (f) intervals of 1 arcmin, across the 
equator. The temperature field is obtained using SYNFAST (see Gorski et al, 2005] ) 



with NSIDE = 4096 and attributing to the grid nodes the temperature value of the 
nearest HEALPIX map cell center. No telescope lobe nor measurement noise was intro- 
duced. The separation angle is computed only from the 9 and <ft differences which allows 
fast computation of the angular correlation histogram. We have checked that, for the 
small maps under consideration, this procedure does not introduce any bias by simulating 
maps with the same area but largely elongated along the equator. Reconstructed angular 
correlation spectra are shown in figure HH] as functions of 7, not cos(^), because the high 
angular momentum information is at small values of 7. The spectra are very different, 
one from the other, due to the large low I value variance of the G\. 

Reconstructing the C\ spectrum from these correlation histograms requires some care. 
First, at high I values, the Legendre polynomials oscillate rapidly and the calculation of 
the integral ( IV. 291) must be performed with a very small angular step. For that pur- 
pose, the correlation spectrum histograms have been oversampled by a factor 10 and the 
intermediate correlation function values obtained from a third order spline calculation. 
Secondly, the angular step of the correlation histogram must be smaller than the map 
step, otherwise some information is lost. For instance, at small separation angle, it is 
necessary to distinguish the angular separation angle of adjacent cells from the one of 
cells which are neighbours on a diagonal. We have used a sampling four times finer than 
the interval between adjacent map cells (\ arcmin). 

With these cautions, reconstructed Ci are shown in figure fT71 and fT8l The filter shapes 
have been chosen according to the discussion of the preceding section. Figure [UJ shows 
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the average of the reconstructed Ci spectra over the 173 simulated maps, while figure [181 
shows one particular simulation. Figure [T7] compares also the calculated variance of the 
Ci (formula IV.27f) with the variance deduced from the 173 simulations. They agree well 
for I < 1500. The discrepancy at larger I values may be due to the HEALPIX map 
resolution, and to the histogram binwidth definition. This is reflected by the systematic 
bias of the average reconstructed C\ spectrum for I > 2200. Figure [18] compares the 
reconstructed C/ spectra obtained with two different apodisation cut-off parameters. The 
one corresponding to 9 = 5 deg looks better than the one obtained with 9 = 10 deg, but 
this reflects only the fact that the C\ values are correlated on a larger scale in the first 
case, washing out some fluctuations, as shown in figure [191 
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Figure 13: Predicted < Cf > from £(7) integration for a spherical cap of aperture 

- the red curve shows the predicted 1(1 + 1) < Cf > /2n 

- the black curve the input 1(1 + l)Q/27r spectrum. 
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Figure 14: Predicted 1(1 + 1) < Cf > from £(7) integration with apodization for a 
spherical cap of aperture 19 deg: 

• Top figure 

- the yellow curve is for apodization (9q = 15 deg, A = 0.5 deg) 

- the blue curve is for apodization (6> = 10 deg, A = 5 deg) 

- the black curve is for apodization (9 = 1 deg, A = 2 deg) 

- the red curve is the 1(1 + 1)Ci/2tt input spectrum (the curve is undistinghishable from 
the yellow one at large /). 

• Bottom figure 

- the black curve is for apodization (6*0 = 10 deg, A = 2 deg) 

- the red curve is the 1(1+ l)Ci/2ir input spectrum (the curve is not plotted above I = 200 
because it is undistinghishable from the black one). 
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Figure 15: Predicted variance and correlation for Cf with apodization for a spherical 
cap of aperture 19 deg: 

• Top figure: 1(1 + l)ergf/27r 

- the black curve is for apodization (6*0 = 10 deg, A = 2 deg) 

- the red curve is for apodization (9$ = 10 deg, A = 5 deg) 

• Bottom figure: correlation values for I = 1500 

- the black curve is for apodization (9q = 10 deg, A = 2 deg) 

- the blue curve is for apodization (9 = 5 deg, A = 2 deg) 
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Figure 16: Reconstructed angular correlation function from 173 simulations of 17 deg x 
17 deg squared maps: 

• Left figure 

- the yellow, orange, green and red curves are reconstructed angular function on 4 indi- 
vidual generations. 

- the black curve is the mean angular function computed from the 173 generations. 

- the blue curve is the dispersion of the angular functions computed from the 173 gener- 
ations. 

• Right figure 

- the black curve is the mean angular function computed from the 173 generations. 

- the red curve is the expected angular function computed directly from the input C\ 
spectrum 
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Figure 17: Reconstructed 1(1 + 1) < Cf > /2n from £(7) integration with apodization on 
173 simulations of 17 deg x 17 deg squared maps: 

- the black curve shows the mean reconstructed 1(1 + 1) < Cf > /2n for apodization 
(0 O = 10 deg, A = 2 deg) 

- the red curve shows the 1(1 + 1)Ci/2tt input spectrum 

- the blue curve shows the dispersion of the reconstructed Cf on the generations for 
the same apodization 

- the orange curve shows the dispersion cr^ predicted by ( IV. 271) and ( 1V.30I) for a 17 x 
17 deg 2 spherical cap 
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Figure 18: Reconstructed l(l + l)Cf /2tt for integration with apodization for one individual 
generation of a 17 deg x 17 deg squared map: 

- the red curve is for apodization (9q = 10 deg, A = 2 deg) 

- the blue curve is for apodization (9q = 5 deg, A = 2 deg) 

- the black curve shows the 1(1 + 1)Ci/2tc input spectrum 
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Figure 19: Reconstructed correlation of Cf for integration with apodization for L = 1500: 

- Top figure is for apodization (6q = 10 deg, A = 2 cfeg) 

- Bottom figure is for apodization (#o = 5 deg, A = 2 deg) 

- blacks curves are the reconstructed Cor(l, L) on the 173 generations of 17 deg x 17 deg 
squared maps 

- red curves are the predicted Cor(l,L) for a 17 x 17 deg 2 spherical cap 
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VI Conclusion. 



In these notes we have computed the biases and covariance of the C\ reconstructed 
from small sky maps. This was done for two methods, the widely used one based on FFT 
analysis of the temperature field, and the one which uses the angular correlation spectrum. 
These two methods are complementary. Estimators of the C\ have been defined and we 
have shown how they can be computed numerically, as well as their covariance matrix 
and correlations. These calculations are simpler to perform in the case of sperical cap 
maps, due to the high degree of symetry. We have shown, using simulated sky maps, 
that most of the results do not depend on the map shape at first order. Most of the 
biases introduced by the use of small maps can be studied numerically without having to 
simulate large amounts of sky maps. We have also studied the complicated dependency of 
the reconstructed Q variance over the map size. In the case of the correlation spectrum 
we have shown how its apodisation works. 
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APPENDICES 



A Computation of the angular correlation integral. 



Let's compute the internal integral X of the angular correlation function of section UT3t 
X= J Y^^Y™ 2 ^) dfi x dfi 2 - cosfr)) 

S 2 xS 2 

We rotate the original coordinate axis Oxyz to OXYZ such that Ox — (8i,<Pi) be the 
unit vector along OZ. Such a rotation is performed by first rotating by an angle 0i 
around Oz. The Oy axis moves to Oy\. Then a rotation with angle Q\ around Oy\ 
is performed. The rotation which transforms Oxyz to OXYZ is the rotation of Euler 

angles R(a — <f>i,(3 = #i,0). Let's (6*2,02) be the coordinates of f2 2 in Oxyz et (02, $2) 
in OXYZ. 

The transformation law of the spherical harmonics is (see Messiah , 1964| ): 



The inverse rotation gives: 



As i? is a unitary transform S^Ji?" 1 ^)) = £^,(#(0,0)) 

so we obtain: 



m' 2 

Let's go back to X. As the transform is a rotation, the jacobian is equal to 1 and: 
dQidfl2 = dcos(9i)d(f>i dcos(9 2 )d(j)2 = dcos(9i)d(f>i dcos(02)<i$2 

The scalar product is invariant under rotation, so the argument of the delta distribution 
remains unchanged: Qi.Q 2 = cos(©2). 



S 2 xS 2 m' 2 2 

5(cos(92) — cos(7)) 



= / itc^^je C(7^ 2 ) ©suw^^o) ^os(^)#i *h 

S 2 xS 2 mi 
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Since ^( 7 ,* 2 ) = ^SSI^H^e'*, 

for the integral over $ 2 to be non zero, we must have m' 2 = 0. 



1= J Y^ie.Ai) ^(oos( 7 )) ©aW 1 ,fl 1 ))27rdoo8(W 1 

with 

and P;°(cos(7)) = P/ 2 (cos(7)), we obtain: 

1 = J YFtfufa) ^ l Y^{e 1 ,<t> 1 )J^P l2 {cos{ 1 ))2'Kdcos{e 1 )d<t> 1 
= 2tt 5 hh 5 mim2 P 2 (cos(7)) 



mirri2 P 2 (cos( 7 )) (A.l) 
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B Computation of the C\ covariance for a portion 
of sphere. 

We have seen (cf section liii.2l) that we can write: 

a>im — Bi m .i> m i a>i> m i 

I'm' 

with either 

um ^r h / (2z + i)(2f' + i)(2f» + i) /z' z r\/z' z z«\ 



i",m" 

or 



Blm^l'm' — J Y™ Y™ dVt 



A 



The relation Yf m = (-l) m Y™ leads to: 

Blm-l'm' = Bi'm'-jm = ( — l) m+m S / _ m;i /„ m / £ C 

Let's compute the ensemble average < Cj >: 

^ 2Z+1 5E ^Zm C^Zm 

m=— I 

= 2I+T S E E Bl m -ir m r Blm-l'/m'i a\i m i ai" m " 
m I'm' l"m" 



< Ci > — 2I+T E E E Blm;l'm' Blmp'm" < a /'m' Ct/"m" > 
m Z'm' i"m" 

7T7 E E Yl Blm;l'm' Bl m yi m n 5l'l"Sm'm" Cl 



21+1 

m I'm' l"m" 



<Cl> ~ 2/~+T 2 I 



+1 

2 a 



m=—l I'm' 



Now let's compute the C\ covariance: 

Ci Cl = (m+i)(2l+i) EE E E E E 



m Z'm' Z"m" M L'Af ' L" M" 



Blm;l'm'Bl m -l" m " B LM;L> 'M' B 'lM;L" M" Q>l'm'0>l"m" &L' M' M" 

Using the ensemble average of 4 a/ m products for a real gaussian temperature field: 

^ ^Ziroi^^m^ ^£3^13 ^47714 ^l\^lz ^Ixli^mxm^l^li^m^m^ &l\l/fim\m$l?X^m2mo, t 

_|_ J_|jmi+m2 G\ x Ci % 6i 1 i 3 6 mi - m3 6i 2 i 4 6 m2 - m4 
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We obtain: 

< Ci Cl > = (2Z+l)(2L+l) E ' ' ' E Blm;l'm'Bi m; i»m" B LM ;L> M> B LM -L"M" 

{Cl>C L i (Syy>S m i m '/S L ' L "S M 'M" + $l>L"$m'M"fil"L'$m"M>) 
+ ( — l) m+m CyCyi Sl'L'Sm'-M'8l"L"8m"-M"} 

<qc l > = {2l+1) \ 2L+1) EE{ 

m M 

E CyCjj Bi m -y m >Bi m -y m i B LM;L' M' B 'lM;L> M> 
I'm'L'M' 

+ E Cl'Cv Bi m .y m >Bi m . L i M i B LM . L i M iB LM .y m i 
I'm'L'M' 

+ ( — l) m+m Yl CyCyi Bi m .y m 'Bi m .y/ m ii B LM .y_ m iB LM .yi_ m '/ 
I'm'V'm" 

} 

The first term corresponds to < C\ >< Cl >■ 

One could exchange l",m" to V , M' in the third term sum: 

V(l,L) = <C l ~C~ L > - <C l ><C L > 

= (22+l)(2L+l) SEE E CyC V { 
m M I'm' L'M' 

Blm-J'm'Blm-L'M' B LM;L> M' BlM-I'tu' 

_|_ ^_^m+M Bi m -jim'Bl m -L'M' B LM;l> -m> B LM;L' -M> 
} 

Then exchanging M to — M in the second term sum and using the symetry relation for 
the Bi m .y m r. 

V ( l i L ) = (22+l)(2L+l) E E E E C v Gu { 

m M I'm' L'M' 

Blm-l'm'Blm-L'M' B LM;L' M> B LM;l'm' 

_|_ ^_^m+M Bi m -y m /Bi m -L'M' B L-M;l'-m' Bl-M;L'-M' 
} 

= (2j+i)(2L+i) E E E E CvCv { 

m M I'm' L'M' 

Blm-l'm' B l m -L'M' B L M;L' M > B LM .y m i 
+ Bi m .y m >Bi m . L i M > B LM .y m 'B LM . L ' M > 

} 

= (imfeny E E E E c i ,c v { 

m M I'm' L'M' 

Blm-l'm' Bl m -L'M' B L M;L'M' BlM-I'tu' 
} 

— (2l+l) 2 (2L+l) EE (EC/' Bi m .y m iB L M;l'm' ) I E C L i B lm -U M' B L M;L'M' ) 
m M \l'm' J \L'M' J_ 

= (2l+l) 2 (2L+l) EE (EC/' Blm-l'm' BLM;l'rn' I ( E ^V Bi m -L'M' B L M;L'M' I 
m M \l'm' J \L'M' J 

We finaly obtain: 

V(l,L) = <CiC~ L > - <Q><Cl> 

+1 +L 

= (2l+l)(2L+l) E E I EC' Blm-l'm' B L M;l'rn> \ 
m=-l M=-L I'm' 
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If we replace the Bi m .i> m > by their values, we obtain: 



V(l,L) 



(2/+l)(2L+l) 



+1 +L 

E E 

m=— I M=— . 



E c v 

I'm' 



X 



X 



,M 



(47T) : 



E hp (-1) 

6/9 
.d«J 

E E 

m=-J M=-L 

V 

m' 



(2f+l)(2F+l)(26+l) 

47T 



(2L+l)(2F+l)(2d+l) 

47T 



to' 
m' 



Z b 

-m (3 

L d 

-M 5 



V I b 


V L d 




E (2/' + 1) C v x EE y/W+Wd + T) bbpbas 

I'm' b/3 dS 



X 



I b 

—m (3 



I' 

m! 



L d 
-M 5 



V I b 




V L d 




That formula is not very usefull because it involves the computation of a general 3j symbol 
which is time consuming. This is why we have to add another constraint: the spherical 
symetry of the portion of sphere under study. 

Note that we would have obtain exactly the same result by doing the computation with 
the other definition for the Bi m y m i. In this case, because it would involve integral of 
product of Y™ other a portion of sphere of general shape, the computation time would 
be enormous . 
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C Computation of the angular correlation distribution 



for a portion of sphere. 

Let's compute 
N A (l) £(7) = E Oh^b LlMl a hm2 b L2M2 

S 2 xS 2 

where the sum ^ runs on mi, Li, Mi, £2, "7,2, -^2, M 2 }. 

We perform the integration on Q 2 by rotating the frame (Oxyz) into (Ox'y'z') such that 
the axis Oz' be on fli (see figure [20]) 




> 

y 



Figure 20: Euler rotation. 



We perform the Euler rotation R(a,/3, 7) of the framed): 

• Rotation of a = (pi around Oz: (0,x,y,z) — > (0,xi,yi, Z\ — z) 

• Rotation of (3 = 9 1 around Oy\. (O, Xi, yx, Z\ — z) — > (O, x', y' = yi, z') 

• Rotation of 7 = around Ozi (i.e. the identity). 
So the rotation is: R(a, /3,j) = R((f>i, #i, 0). 

We define (8, 4>) to be the polar coordinates of the frame (0, x, y, z) and (8', </>') the polar 

6 Here 7 is the third angle of the Euler rotation not the separation angle of the angular correlation 
function. 
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coordinates in the new frame (0,x',y',z'). 



,(0 

ra'ra v 
m! 

mm' 



= ESJL^iA.o)!?"'^) 

Performing the rotation we get: 



A/" A (7H(7) = E^h lMl a hm2 b L2M2 f d^Y^i^Y^in,) 

S 2 



xj dn> 5(n[.Q' 2 - oo8( 7 )) s^C' 2 (^)®S 2 M^Lf(^ 2 ) 

where the sum runs on {li, mi, Li, Mi, l 2 , m 2 , L 2 , M 2 , m 2 , M 2 } 
and where we simplify the notation ®^, m = S3^/ m (0i, #1, 0). 

We have dfi' 2 = dcos(9' 2 )d(f) 2 . The integral on dcos(9' 2 ) and the function £(•••) leads to 
the replacement of 9' 2 by 7. Writing Y l m (9,(f)) = Xf l (9)e im ' t ', we have for the integral on 
'■>'/■ 

f J d cos( W 2 (fi 2 )F L f (Q' 2 ) = J d<\>' 2 \£ (7) \f 2 (7) e^+^ 

= 27r<5 m ,,_ M , a£( 7 )aJJ( 7 ) 

Thus we obtain: 

^(7)^(7) = 2vrE a: 2 (7)az: 2 (7) 



5 2 



where the sum runs on {li, mi, Li, Mi, l 2 , m 2 , L 2 , M 2 , m 2 }. 

Now we perform the ensemble average, and as < &i imi cii 2rn2 > — C^S^^S^^, we get: 



< A^( 7 K(7) > = 2n Z b LlMl b L2 M 2 C h A™ 2 (7^ (7) 



5 2 



where the sum runs on {li, mi, Li, Mi, L 2 , M 2 , m 2 }. We have: 



and by definition, ft; is on Oz', (0[ = 0), so Y™ 2 ^) = A"^ 2 (0)e" im ^i . 
As A^(0) = A° (0)(J^o, we obtain: 

< A/^( 7 ) £(7) > = 2tt E W>l 2 m 2 C h A°(7)A° 2 ( 7 )A°(0) J dtotfft&fa, 6 U 0)Y^(^i) 

s 2 

52 



where the sum runs on {li, Li, M 1; L 2 , M 2 }. 



But we have: ft, 0) = J^Ygfak. 



<N A { 1 )£{ 1 )> = 2nJ2J J ^ ri b LlMl b L2 M 2 C h \l( 1 )Xl( 1 )Xl(0) 



xj dQ 1 Y^(9 1 ,cf> 1 )Y^(Q 1 ) 
s 2 



2vr £ JjSti h lMl b L2 M 2 C h A°(7)A° 2 ( 7 )A°(0) 8 L2Ll 8 AhMl 



= 2vr £ ^/hJtt WAxM, C h A? i (7)A Li ( 7 )A°(0) 
where the final sum runs on Li, Mi}. 



Using the definition of B x (see [XTLMj) as well as \f = ^^W|j=^Ff», P° = Pi and 
Pi(0) = 1, we obtain: 

< M A iri) £(7) > = (e ^ c h p h (^ x (2vr e (2L t + i) b Li p Li ( 7 )) 



To compute 7V A ( 7 ) = / dfi lC /fi 2 W A (^i) W A (fi 2 ) 5(^1-^2 - cos( 7 )) we do the 

5 2 x5 2 

same computation as before with T(Q) = 1 = V^7rF °(fi) so ai imi = V^br^o^miO and 
C h = Att8i i0 ■ 

We thus obtain (P (6) = 1): 

M A {i) = 2nJ2 (2£ + 1) Pl(7) 

and finally: 



<^)>= E^^W 



2l + l^^A^ i2L+1)BLPM 



471 J E ( 2 -^ + 1) b£ Pl{i) 

L 



53 



D Computation of the angular correlation integral 



for a portion of sphere. 

Let's compute the integral: 
i*£,(c) = / dVtdVt' 5{Q.Q' - c) Y^' (Q')Y™(tt) 

AxA 

= J dttdtt' Sitt.tt' - c) W A (9!)W A (n) Yf (tf)Y{n(n) 

S 2 xS 2 

E bLMh'M' J dQdQ' 5{nS\' -c)Yp\Q, , )Ylf\n ! )Yf\n)Y^(Sl) 

LML'M' S 2 xS 2 

Using the formula of the product of two Y™ relative to the Clebsch-Gordan coefficients 
and remembering that the coefficients are real: 



Jim ( r \ _ ST h h V / 2i+l 2L+1 / (2f'+l 2L'+1 

WW ~ LM b L > M > 1^ A/ 4 ^(2A+1) V 47r(2A'+l) 

x < ILmM | Xfj, >< /MOO | AO > < I'L'm'M' | A'// >< Z'L'OO | A'O > 

x / dQdQ' 5(Q.Q' - c) Y$'(tt)Y£(n) 
s 2 

The remaining integral has been computed in appendix [X] 

X = 2tt 5 x\' <W P\(c) 

Thus 



l\Z'{c) = 2tt E b LM b UM > E ^a(c) 

LML'M' Xfi 



x W ^sa+i^^ < I V >< iMOO | AO > 

(2i'+l)(2L'+l)' ^ VTI^I^/Tl 



x V M2A+1) < l L m M | A/i >< Z'L'OO | AO > 
It is real and can be rewritten as integrals of three Y™. 



LML'M' A^ 52 5 2 

= 2ttEPa(c) E /dfiif r (n)&^(n)^(n)/dn'ir'("OWV^'(nO>T(o') 

A/i LML'M' s 2 S 2 

= 2ttEPa(c) / rffi ^(f2)y L M (n)y^(n) / do' w^(tf)Y#'(tf)Y?(fi') 

Am S 2 S 2 

= 27rE^A(c) / dQYf{Q)Y x u {Q) j dD! Y™' (Vt')Yf(Q!) 

Xfi AxA AxA 

Remembering the Bi m ii m i definitions (see section liii.2l) : 



LMl'm' 

LM 
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